arXiv:1508.04259v4 [cond-mat.mes-hall] 5 Jan 2017 


Charge-insensitive single-atom spin-orbit qubit in silicon 

Joe Salfi/’^ Jan A. Dimitrie Culcer/ and Sven Rogge^’^ 

^School of Physics, The University of New South Wales, Sydney, NSW 2052, Australia. 
^Centre for Quantum Computation and Communication Technology, 

The University of New South Wales, Sydney, NSW 2052, Australia. 

(Dated: January 6, 2017) 

High fidelity entanglement of an on-chip array of spin qubits poses many challenges. Spin-orbit 
coupling (SOC) can ease some of these challenges by enabling long-ranged entanglement via electric 
dipole-dipole interactions, microwave photons, or phonons. However, SOC exposes conventional spin 
qubits to decoherence from electrical noise. Here we propose an acceptor-based spin-orbit qubit in 
silicon offering long-range entanglement at a sweet spot where the qubit is protected from electrical 
noise. The qubit relies on quadrupolar SOC with the interface and gate potentials. As required 
for surface codes, 10® electrically mediated single-qubit and 10^ dipole-dipole mediated two-qubit 
gates are possible in the predicted spin lifetime. Moreover, circuit quantum electrodynamics with 
single spins is feasible, including dispersive readout, cavity-mediated entanglement, and spin-photon 
entanglement. An industrially relevant silicon-based platform is employed. 

PACS numbers: 71.70.Ej,73.21.La,42.50.Pq,03.67.-a,03.67.Lx 


In recent years, the coherence and control fidelity of 
solid-state qubits has dramatically improved [TH5] and 
spin qubits [SHE] with highly desirable properties have 
been demonstrated. lanni However, many obstacles re¬ 
main to efficiently entangle a large array of spin qubits on 
a chip. For example, exchange is inherently vulnerable 
to decoherence from electrical fluctuations [1 lIfT5] . cou¬ 
pling spin to charge noise. Minimizing decoherence and 
improving control in the face of noise is the key issue for 
large-scale quantum computing, because it ultimately de¬ 
termines if the error-correction resources can be managed 
for a large qubit array. m Moreover, exchange-based en¬ 
tanglement is inherently short-ranged, making fabrica¬ 
tion challenging for gates in quantum dot arrays [5], and 
placing strict demands on Si:P donor placement. [7] 

Here we propose a single-acceptor spin-orbit qubit 
where the unique properties of hole spins give a host 
of desirable attributes. First, spin-orbit coupling (SOC) 
enables long-ranged entanglement via microwave pho¬ 
tons or electric dipole-dipole interactions . of inter¬ 
est for hybrid quantum systems . improving error 

correction]^, and reducing fabrication demands com¬ 
pared with exchange coupled schemes. Second, and 
most remarkably, we find a sweet spot where coherence 
is insensitive to electrical noise and electric dipole spin 
resonance [35H33 (EDSR) is maximized. Consequently, 
coherence and gate timings are protected from electri¬ 
cal noise at the Hamiltonian level, and one- and two- 
qubit gate times are optimized. In comparison, electric 
field noise dephases conventional spin-orbit qubits [551? 
] and acceptor charge qubits. Ha [33 The coherence of 
our spin-orbit qubit benefits from reduced hyperflne cou¬ 
pling of holes [3S] and ^®Si enrichment [33] , and has much 
longer phonon relaxation times than acceptor charge 
qubits. [231 Ei] Finally, the acceptors naturally confine 
single holes that can be manipulated in silicon nanoelec- 


tronic devices [40]. 

The exceptional properties of the qubit derive from the 
quadrupolar SOC [4Tlj44] contained in the spin-3/2 Lut- 
tinger Hamiltonian[45| and in the interaction with the 
inversion asymmetric interface potential, not studied pre¬ 
viously for acceptors. This SOC is unusually strong for 
acceptors because it acts directly on the low-energy spin 
manifold, contrasting its indirect role in hole quantum 
dots.[T9l [20l I46ll49] The SOC must be considered non- 
perturbatively to obtain the sweet spot, and the interface 
strongly enhances EDSR relative to a bulk acceptor. We 
find 0.2 ns one-qubit gate times, charge-noise immunity, 
and long phonon relaxation times at the sweet spot, al¬ 
lowing for > 10® operations in the coherence time. Two- 
qubit entanglement based on spin-dependent electric 
dipole-dipole interactions [T5HT7] is feasible with -s/SWAP 
times of 2 ns, and 10^ operations in the coherence time. 
EDSR also enables circuit quantum electrodynamics [551- 
l30] (cQED) with single-spin dispersive readout, and long 
distance spin-spin entanglement with -s/SWAP times of 
200 ns. Resonant spin-photon coupling with = 5 MHz 
is also feasible. 

Qubit Concept. The qubit is a hole spin bound to a 
single Si:B dopant [40[ [50l [51] ; implanted [52] or placed by 
scanning tunneling microscopy |53l [51] near an interface, 
in a strained silicon-on-insulator (SOI) substrate (Fig. 
m The key quadrupolar interactions, associated with 
interface inversion asymmetry and products {Ji, Jj} = 
(JiJj + JjJi) of spin-3/2 matrices where i[j) = x,y,z, 
originate from strong SOC in the valence band, and have 
no analog in the conduction band. [lTH5i] This SOC acts 
on the 4x4 ground state manifold be., the mj = 

±1 and mj = ±5 Kramers doublets composed mostly 
of |j = ^,mj) Bloch states. [55] For Si:B they are well 
isolated by ~ 20 meV from orbital excited states and 46 
meV from the valence band edge[35] (Fig. m- 
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The key quadrupolar interactions include the accep¬ 
tor hole spin-mixing that is linear in electric fields, 
^^E,ion = ‘^p/Jy} + c.p), associated with Td 
symmetry in the central cell [57]. Here, p = 0.26 D is 
known for SirBjSSj (1 D = 0.021 e-nm). An electric field 
Ez further breaks the envelope function parity by mix¬ 
ing excited states outside the manifold. |55| Pro¬ 
jected into the subspace, this interaction is gov¬ 

erned by He = b{j‘^ - \l)El + {2dlVY>){{Jy, Jz}EyEz + 
{Jz, Jx}EzEx), where b and d split and mix the dou¬ 
blets, respectively. We verified that this holds for 
triangular interface wells, using (i) a Schrieffer-Wolff 
transformation |59l [60] with higher excited states in the 
spherical spin-3/2 basis[61). and (ii) numerical, non- 
perturbative Luttinger-Kohn (LK) calculations with ex¬ 
plicit ion and interface well potentials [5U [53]. We find 
a splitting AwiEz) = Ait + Ag^Ez) (Fig. [^), where 
Ait from the interface is larger for shallower acceptors 
(in agreement with experiments |50j ). and Ag{Ez) oc Ez 
increases with increasing field. Moreover, quadrupolar 
SOC combining inversion asymmetry and in-plane elec¬ 
tric fields is governed by terms a{Ez)Ex,y oc EzE^^y that 
replace dEzE^^y hi He- 

Operating point and sweet spot. Here we show that the 
qubit splitting huj (between |-|-) and |—), Fig. IB) in an 
in-plane magnetic field yB depends on the electric field 
Ez applied by the gates (Fig. lA), and at the sweet spot, 
hoj is insensitive to electric-field noise dE in all direc¬ 
tions. Including magnetic fields, strain A^ (Supplemen¬ 
tal Material]!!]) and the interface well, but not in-plane 
electric fields, we find an operating point Hamiltonian, 

/ A{Ez) -iez i^£z -ipEz \ 

^ ^ i£z ^{Ez) ipEz -i^£z qx 

V ipEz i^£z 0 0 / 
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FIG. 1. A. Device schematic, showing near-interface Si:B im¬ 
purity with gates to SG and TG to apply in-plane and vertical 
electric fields, respectively (left), or for cQED, gates form¬ 
ing the resonator apply both the in-plane and vertical elec¬ 
tric fields (right). An in-plane applied magnetic field ensures 
a long photon lifetime in the superconductor resonator. B. 
Electronic structure of an acceptor, where the splitting A is 
determined by the strain, interface, and gate field Ez . Shown: 
pEz-induced mixing of states in the 4x4 manifold due to the 
Td symmetry in the unit cell of the ion. Not shown: LH-HH 
coupling from in-plane drive fields. 


quencies. We treat the quadrupolar SOC term pEz 
by a rotation that maps pEz exactly to the diago¬ 
nal, to a basis {|^—), \l +), |u—), |u-|-)} leaving Zeeman 
terms off-diagonal. We obtain |1±) = ol \^±i/ 2 ) ± 
ioH |'I';p 3 / 2 ), a low-energy Kramers pair (energy £/ = 
I (A — •\/A2~T4 E|p^)), and an excited Kramers pair 
|n ±) = ql |^± 3 / 2 ) T ian I^Ti/ 2 ) (ener gy £„ = | (A -h 
a/A^ -h AE^p'^)). Here, op = £i/y/E^p^ + ef and 
Oh = a/I — a| = Ezp! E^p^ In the basis 

{\l—) , |/-l-), |m—) , |u-|-)} Eq. 1 becomes 


in the basis {|^'-i/ 2 ) , |^'i/ 2 ) , |^- 3 / 2 ) , |^ 3 / 2 )}, where 
£2 = gip-sB, pb is the Bohr magneton, gi = 1.07 is the 
Lande g-factor for Si:B.[5S], and A{Ez) = Aw{Ez) — A^ 
is the splitting between the light and heavy holes. The 
cubic g-factor [58] 32 ^ ffi is temporarily neglected. 

Inspecting H^p, Ez mixes |'I'±i/ 2 ) and |'I'=p 3 / 2 ) and 
these states have an avoided crossing when the interface 
well splitting compensates strain, i.e., A(if°) = 0. In 
Fig. lA we show that for appropriate strains A^ > Ajf, 
the anti-crossing can be obtained at E^. ^ 15 MV/m for 
zo ~ 5 nm acceptor depths. 

The field E^. at such an anti-crossing is large enough 
that the level-repulsion gap Agap = exceeds the 

Zeeman interactions, i.e., £z/Agap ~ 0.1. This un¬ 
usual aspect of our hole spin-orbit qubit c.f. other 
proposals [T5l - f20| follows from the tunability of the spin- 
3/2 levels with strain and confinement, giving rise 
to the anti-crossing, and the strength of quadrupo¬ 
lar SOC [55] relative to typical spin qubit Larmor fre- 
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Here, the Zeeman terms Xzi depend explicitly on Ez due 
to the gate-induced mixing of |'I'±i/ 2 ) and |4' 3 = 3 / 2 ). We 
find Xzi = 2eziV^aLaH - ia\), Xzu = 2eziV^aLaH - 
iajj) and Xzo = 2£z{-aHaL + i^a\ - i^ajj). 

We perform a final rotation that exactly maps Xzi 
and Xzu to the diagonal, leaving Xzo off-diagonal, 
defining a basis {|—), |+), |e—), |e-|-)} (see Supplemen¬ 
tal Material]!!]). To zeroth order in XzojiTu — £ 1 ), the 
splitting of the Kramers pair qubit states |-|-) and |—) 
is huj « \Xzi\- When mixed by the gate electric field, 
the spin 1/2 and spin 3/2 states with different Zeeman 
terms define a qubit |±) where hio is maximized (inde¬ 
pendent of electric fluctuations) zdEz to first order when 
|Z±) = ^ l^'ii/a) ± ^(- 5 ) |^'=p 3 / 2 ) (see Supplemental 
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FIG. 2. Spin qubit levels £± and eu± for (A) zq = 4.6 nm and (B) zo = 6.9 nm, to zeroth order in Xzoji^u. — Si)- Qubit 
frequency for (C) zo = 4.6 nm and (D) zq = 6.9 nm using approximate (black), analytic (green), and full numerical (blue 
squares) models, in Bq = 0.5 T. Spectral weights |a/,| (blue dashed) and |aj^| (red dashed) are shown. EDSR coupling D for 
(E) Zo = 4.6 nm and (F) zo = 6.9 nm. We take Ae = 0.62 meV (0.34 meV) for zq = 4.6 nm (6.9 nm) achievable in SOI|65], and 
exceeding disorder strain [^ I66j . Parameters Aif, Ag{Ez), and a{Ez) were obtained non-perturbatively in a 6 x 6 LK basis 
including the cubic LK terms and the split-off holes. 


Material [M]). As we will subsequently show, the qubit 
is also insensitive to in-plane electric noise SE^^y, while 
a similar analysis yields another sweet spot at Ez = 0- 
Energy levels e± = ei± for the qubit are shown 

alongside excited levels £e± = ± for zq = 4.6 

nm (6.9 nm) in Fig. [2]A (Fig. [^). Here, blue (red) hue 
denotes the amplitude of ql (an)- The qubit frequency 
is shown in Fig. HP and Fig. HP for approximate (black) 
and exact (green) solutions to iLop, alongside the numer¬ 
ics (squares). The maxima in huj in Fig. HP (Fig- HP) 
defines the sweet spot at Ez = 17 MV/m (14.8 MV/m), 
for loLp = 3/4, as expected. We note that the approx¬ 
imate solution (Fig. HP ,D, black lines) captures the es¬ 
sential behaviour of the analytic model (Fig.[^,D, green 
lines). Corrections to Zeeman interactions from interface 
inversion asymmetry and cubic Lande g-factor, although 
included in the numerics (squares), have been neglected 
in the analytic model (green). Note that the interface 
prevents ionization; although Ez ^ 15 MV/m is much 
smaller than silicon’s breakdown field, it well exceeds the 
ionization field of Si:B. m 

In-plane electric fields: EDSR and noise immunity. We 
express interactions with in-plane electric fields in the 
basis {|-), |-b) , |e-), |e-b)}, yielding 

( Si — ^ 0 oiEi \zi oiEo -b Az2 \ 

- 0 Si E ^ oF/ 2 -t- Az2 oiE\ -b \zx 

aEl+\*z^ aEi + \l^ q 

\ -b \*z^ aEi -b Xzi 0 Eu + / 

(3) 

Here, |-b) and |—) are our Kramers pair qubit states. 


Xzi oc Xzo and A^^ ot Xzo are Zeeman terms, and £’ 1,2 
are interaction terms with in-plane electric fields, where 
El = i(sin0 -b 7] cos 9)Ex + i{cos0 -b r]siii9)Ey, £2 = 
(—COS0 -b r]sm9)Ex + (sin0 — r]cos9)Ey, 9 = 9^ — 9i, 
>^Zi = \Ezi\ exp(z 6 »i), and ry = p/a. 

The qubit Hamiltonian iLqbt = Ejjaz + ££||CT 2 ;, where 
hu! is the qubit frequency (Fig. 2C,D) and D is the EDSR 
matrix element (Fig. 2E,F), is obtained by projecting the 
off-diagonal elements of H to first order in Ex,y using 
a Schrieffer-Wolff transformation. [59l l60] Notably, qubit 
coherence is protected from in-plane electric noise since 
hoj contains no terms to first order in Ej;y. EDSR drive 
comes from the transverse coupling DE^^a^ in i7qbt- We 
obtain D = a\Xzo\{£i — £u)~^{o cos{9o — 0 ||) -bpsin(0o — 
0||)), where Ey = £||(xcos0|| -bysin^y). Interestingly, 
the small splitting £„ — £/ essential for spin mixing at the 
sweet spot also causes strong EDSR, since D oc (£„ — 
£i)~^. Note that the EDSR term is dominated by the 
inversion asymmetry quadrupolar SOC parameter a « 
25 D (Fig.[2p), since it is 100 x larger than the bare Td 
SOC parameter p. 

Importantly, D can be maximized at the sweet spot by 
choosing the angle Ey relative to B||y (see Fig. 2E,F). 
This yields fast gate times, but it also makes D, and 
therefore all timings based on EDSR, insensitive to fluc¬ 
tuations in electric field, protecting gate fidelity from 
noise at the Hamiltonian level. Since p = p/a ^ 0.01 
and 9o = 7 r /4 at the sweet spot, D is maximized with 
respect to 0y at 0y = — 7 r/ 4 ± 7 r/ 2 . As shown for zq = 4.6 
nm (6.9 nm) in Fig. (Fig. HP) D is maximized with 
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respect to for the same choice 9\\. This result can be 
easily obtained analytically, and holds for the analytic 
(green) and numerical (blue squares) solutions. 

Qubit Operation. The one-qubit and two-qubit gates 
employ EDSR-mediated interactions at the sweet spot, 
where coherence is protected from noise, and their times 
T are minimized and also insensitive to electrical noise. 
EDSR driven tt rotations require ti = h/{2DEAc) = 1 
ns (0.2 ns) for the zq = 4.6 nm (6.9 nm) deep acceptor, 
assuming a modest in-plane microwave field E^c = 500 
V/cm. A 7 r /2 (0) phase shift realizes a cjy {ax) gate, 
and Gz gates can be decomposed into a sequence of Gx 
and Gy gates. Readout can be accomplished by energy- 
dependent [68] or spin-dependent [69| tunneling, or disper¬ 
sive readout in cQED.|55| Initialization can be achieved 
by projective readout followed by spin rotation. 

Two-qubit entanglement can be achieved via long- 
ranged Coulomb interactions, owing to spin-dependent 
electric dipole-dipole interactions. [IMTT] Their strength 
is given by Jdd = (vi • V 2 i?^ - 3(vi • R)(v 2 • R))/ 47 rei?'^, 
where R is the inter-qubit displacement and is a 
spin-dependent charge dipole of qubit i, which has the 
same magnitude as the EDSR matrix element. Eor a 
20 nm distance with negligible tunnel coupling, we ob¬ 
tain a Vs WAP time of T 2 dd = h/AJdd « 2 ns with 
Jdd ~ /AneR^. The 10^ times enhancement of EDSR 
from the interface reduces T 2 dd by 10 "^ relative to accep¬ 
tors in bulk silicon, and 10 ® relative to bare magnetic 
dipole-dipole coupling. Entanglement by Heisenberg ex¬ 
change is also possible and exchange is hydrogenic when 
A exceeds J. IS] We note that the advantage that holes 
do not have valley degrees of freedom (TOj which may com¬ 
plicate Heisenberg exchange for electrons in Si. HD 

Circuit QED. Coplanar superconducting microwave 
cavities could be used to implement cQED includ¬ 
ing two-qubit gates, dispersive single-spin readout, and 
strong Jaynes-Cummings coupling on resonance with the 
cavity. We assume a coplanar waveguide res¬ 

onator operating at R = 0.5 T (/ = 15 GHz) and a 
vacuum electric field Eq « 50 V/m. This can be ob¬ 
tained using a tapered resonator gap, or a supercon¬ 
ducting nanowire resonator. m At the sweet spot for 
Zq = 4.6 nm (6.9 nm), the vacuum Rabi coupling is 
Pc = cDEq = 2 neV (10 neV). 

Eor cavity mediated non-demolition readout and qubit 
coupling, we detune the qubit from the cavity by A = 
45c- [22] Here, the spin state shifts the cavity resonance by 
A/ = g^lA = 0.25 MHz (1.25 MHz) for zq = 4.6 nm (6.9 
nm). The two-qubit VSWAP time is T 2 c = h/AJ^ = 200 
ns for zo = 6.9 nm, determined by the effective spin- 
spin interaction[32] Jc = 2g\lA = 2.5 MHz. Operating 
at zero detuning, spin/photon Rabi oscillations require 
hir/gc = 1 pLS (200 ns). Assuming Q = 10® at Rq = 0-5 T 
in state-of-the-art superconducting cavities [7^ 175] PcH = 
6.7 (33) Rabi cycles can be obtained for zq = 4.6 nm (6.9 
nm), where k = f/Q is the cavity loss rate. 


Relaxation and Dephasing. We consider spin-lattice 
(phonon) relaxation and dephasing from a host of elec¬ 
trical noise sources, and compare them to gate times. 
Since silicon is not piezoelectric, spin-lattice relax¬ 
ation occurs only via the deformation potential [Til 175] . 
Eor temperatures T hLo/ks, the spin relaxation 
time derived in the Supplemental Material [^ follows 
Tf^ = (ha;)®(C'd/ 20 p 7 rh'^)(|Azo|/A) 2 , where |Azo|/A = 
huj/ApEz at the sweet spot, and Cd = 4.9 x 10“^° 
(eV)^(s/m)®. We obtain Ti = 20 ps (5 ps) for zq = 4.6 
nm (6.9 nm) at Rq = 0.5 T that are 100 times longer c.f. 
bulk unstrained silicon at R = 0.5 T. [551117] 

Random fluctuations in qubit splitting h,5uj{t) dephase 
the qubit. The dephasing rate from random telegraph 
signal (RTS) in charge trap occupation is (T^)”^ = 
(i5a;)^Ts/2, where 5hio is qubit frequency shift, and rg is 
the average switching time. [36] We take ts = 10®ti as the 
worst case, since slower fluctuations can be suppressed by 
dynamical decoupling. Assuming a trap 50 nm away, we 
find 5E ^ 2,000 A7/m and a large window of 200,000 
V/m (20, 000 V/m) of gate space where T| > 2Ti at the 
sweet spot, for zq = 4.6 nm (6.9 nm). In comparison, the 
same analysis gives T| ^0.1 ns for acceptor-based charge 
qubits with similar gate times. It is remarkable that in 
comparison, electrical noise has virtually no effect on co¬ 
herence in our spin-orbit qubit, illustrating the advan¬ 
tages of inversion asymmetry and our spin-orbit qubit’s 
sweet spot. We also find that dephasing from Johnson- 
limited gate voltage noise, and from two-level (tunneling) 
systems (TLS), are ~ 10^ and ~ 10^ times weaker, re¬ 
spectively, compared with RTS.[51] There are only a few 
spin resonance experiments on acceptors [551 [7HH5T] . none 
of which feature strain and an interface. [50] We expect 
hyperfine-induced decoherence in "®'*Si to be weak since 
it has only 4.7 % of spin-bearing isotopes and hyper- 
fine interactions are weaker for holes than electrons. [55] 
Meanwhile, ^®Si enrichment could be used to virtually 
eliminate the nuclear bath. [55] 

The insensitivity to Johnson noise and tunneling TLS 
means spin-lattice Ti limits coherence for few (or slow 
enough) traps at Si/Si 02 interfaces. Eor R = 0.5 T, 
ri > 10 "^ single qubit gates, r 2 dd > 10 ® dipole-dipole 
two-qubit gates, and r 2 c ~ 25 cavity-mediated two-qubit 
gates can be achieved in a Ti limited coherence time. 
Therefore while Ti is short compared to donors, many 
gate operations can performed. Since Ti oc a;“®, choosing 
R = 0.25 T increases all ratios favourably to ri > 10®, 
r 2 dd ~ 10^, and r 2 c ~ 50. Since Ti is much longer at 
the Ez = 0 sweet spot, adiabatically sweeping to Rz = 0 
opens a pathway for a long-lived quantum memory. 

Conclusions. The proposed single-acceptor spin-orbit 
qubit exploits the tunability of the J = 3/2 manifold 
of acceptors and the associated quadrupolar SOC arising 
from the ion and interface potential, providing for (i) fast 
one-qubit and long-ranged two-qubit gates (ii) at a sweet 
spot where the qubit phase and all gate timings are in- 
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sensitive to electrical fluctuations, (iii) avoiding entirely 
the need for exchange interactions, (iv) in an industri¬ 
ally relevant silicon platform. 10 ® single-qubit and 10 ^ 
two-qubit gates could be possible in the qubit coherence 
time. Using cQED, dispersive single-spin readout, cavity- 
mediated spin-spin entanglement, and Jaynes-Cummings 
spin-photon entanglement are possible. 

We thank R. Winkler, U. Zuelicke, M. Tong, and T. 
Kobayashi for helpful discussions. JS, JAM and SR ac¬ 
knowledge funding from the ARC Centre of Excellence 
for Quantum Computation and Communication Tech¬ 
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Research Office (W911NF-08-1-0527). DC acknowledges 
funding through the ARC Discovery Project scheme. 

APPENDIX 

INTERACTIONS WITH MAGNETIC FIELDS 
AND STRAIN 

Interactions of acceptor-bound holes with magnetic 
fields and strain are known for acceptor dopants in bulk 
silicon [55l[57] . In the I'l'mj) subspace, interactions with 
magnetic fields B = -|- yBy + zB^ are represented 

by the Hamiltonian 

Hz = fJ-sigiiJxBx + c.p.) -f g 2 {JxBx + c.p)). (AI) 

Here, Ja are J = 3/2 matrices, c.p. refers to cyclic per¬ 
mutations, gi and g 2 are the linear and cubic Lande g- 
factors, and gs is the Bohr magneton. Interactions with 
strain are represented by the Hamiltonian 

=a'Tr[e] -b b'{{JI - \l)exx + c.p.) 

+ (2d7v^)({J,, , Jy\f^xy + C.p. ), (A2) 

where {J^,, Jy} = \{JxJy + JyJx), cl', b' and d' are Bir- 
Pikus deformation potentials [55l (57] . 

ACCEPTOR STATES IN SPHERICAL SPIN-3/2 
BASIS 

In the spherical spin-3/2 basis \L,J;F,mF), where 
F = L -b J is an effective total angular momentum | 6 I|. 
acceptor eigenstates take the form 

\^m.,) = fo{r)\L = 0,J=^-,F=lmF) 

+5o(c) |T = 2, J = |;F = |,tof) , (A3) 

where J = | is an effective total spin, and the spin- 
3/2 spin-orbit interaction has coupled states with AL — 
0, ±2. The /o('c) and goir) are radial envelope wave- 
functions for envelope function spherical harmonics with 
L = 0 and L = 2, respectively. The |L, J;F, tof) are 
found using the Clebsch-Gordan coefficients. 


States outside the 4x4 subspace, given in ref. m, 
contribute to envelope function asymmetry to realize the 
b and d terms in FIe from the main text. The Schrieffer- 
Wolff calculation mixing in these higher excited states 
will be presented in a future paper. 

ANALYTIC MODEL FOR LOW ENERGY 
STATES 

The low-energy holes are described in the basis 

by 
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where 

= gigsB, and 


gi = 1.07 for B:Si[58l [77]. The unitary transform Uq 
that diagonalizes Hop{ez = 0) is 
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in the basis {| 4 (+ 3 / 2 )^+i/ 2 ) , |^-i/ 2 ) , |^(- 3 / 2 )}- 
Here, of = £i/+ ef and an = — a\ = 

HzP/ + ef. Applying this transformation we 

obtain an extended qubit Hamiltonian for the spin qubit 
states and nearest excited spin states 

/ ^e/ ^^*zo ^ 0 \ 

Hop= K „ ^ l\* ■ (A 6 ) 

2'^Zo t) £u 2^Zu 

V 0 \^Zo \^Zu £u J 

in the basis {|Z—) , |/-b), |u—)}, |M-b)}. Here, the effec¬ 
tive Zeeman interactions Xzi appear on the off-diagonal. 
The effective Zeeman interactions depend explicitly on 
the gate electric field due to the pBz-induced mixing of 
|^±i/ 2 ) and |'I';p 3 / 2 ). These off-diagonal Zeeman inter¬ 
actions are 

Xzi =‘^Szi'^SaLaH — ia^), (A7) 

Xzu = 2£z{V^aLaH - iOn), and (A 8 ) 

Xzo = ^ezi-anaL + iV3al/2 - iy/3a%l2). (A9) 

The upper and lower 2x2 blocks of Hop are diagonalized 

by 

/ +e-ie,/2 0 0 \ 

_g+*e ,/2 ^^+^e,/2 0 0 

0 0 _e +*»“/2 y 

(AlO) 


Uzo = 22 
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in the basis {\l—), |^+), \u—), 1 ^+)} where, 


ELECTRIC DIPOLE SPIN RESONANCE 


0 ; = arctan(-\/3a//OL, —a|), and 
9u = arctan(-v/SaLOff, ~o^h)- 


(All) 

(A12) 


The total interaction with in-plane electric helds 
and Ey described hy El e and He, ion is 


Applying this transformation we obtain the following ex¬ 
tended Hamiltonian 

/ £; — 0 Xzi \z 2 \ 

fr _ 0 ei + Xz 2 Xzi 

°p ~ \* c- Azt .1 n 

V ^Z2 ^Z1 0 £« + / 

(A13) 

in the basis {|—), |+), |e—), |e+)}, where, 

Xzi = ^lAzol cos(dz /2 - d „/2 - Oo), (AM) 

Az 2 = 5 |Azo|*sin(d ;/2 - 6»„/2 - do), and (A15) 
9o = arctan(^(a| — a^), anaL)- (A16) 

Our approximate qubit model from the main text takes 
Xzo /(£« — £;) to zeroth order to give qubit states |—) and 
|-|-). In this approximation we obtain a qubit frequency 

huj = \Xzi\ = 2£z\J‘ia\a?fj + a\. (A17) 

A sweet spot occurs when the qubit frequency is insensi¬ 
tive to small fluctuations in electric fields. For an accep¬ 
tor experiencing no static external applied held along x 
and y directions, this occurs for the roots of 

dhw/dEz = 2ez{daL/dEz){i — '^a\)/ij 3 — 2a\. (A18) 

One root (sweet spot) is {aL,aH) = (—^, 5 )- Substi¬ 
tuting ei and A{Ez), we obtain an equivalent condition 

AG{Ez) + Ai{ + 2p/V5Ez = A, (A19) 

for this sweet spot. Another root (sweet spot) occurs at 
the roots of 

daz/dEz = Ey{deildEz){e^ + Eyy^/^, (A20) 

The sweet spot associated with the above root occurs 
at Ez = 0. The lowest order correction to huj due to 
coupling to levels |e±) is easily obtained from 2 "‘^-order 
perturbation theory. 


He,, =a{Ez) 


I 


0 


+P 


0 

V 0 


(^ 

E_ 

0 

0 

\ 

E+ 

0 

0 

0 

0 

0 

0 

-E. 


V 0 

0 

-E+ 

0 

/ 

— iEjz 

0 

0 

\ 

0 


0 

0 

0 


0 

+iEj, 


0 


-iE_ 

0 

/ 


(A22) 


in the basis {|'I'+ 3 / 2 ) , |'I'+i/ 2 ) > |^-i/ 2 ) ; |'I'- 3 / 2 )}j 
where Ej^ = E^ + iEy and E_ = E^ — iEy. The hrst 
matrix is the coupling due to the broken inversion 
symmetry of the interface and gate held, along the 
z direction, while the second matrix describes the 
interaction due to the symmetry of the local held of 
the ion. Rotated into the qubit basis using Uto = UzoUq, 
we obtain 


H = 


0 

aEl + Z* 

\ aE^ + Z 2 aEl + Z* 


aEX + ZX 


aEl 

aE2 

£ 11 . 


“h Z 2 

\^Zu I 


+ ^2 \ 


0 


aE2 Z 2 

aEl + Zi 

0 

I o 


(A23) 

in the basis {|—), |-|-), |e—), |e-|-)}, where, Ei = i(sind-|- 
y cos 9)Ex + i(cosd -I- 'qs\n9)Ey, E 2 = (—cosd -|- 

'qs\\i9)Ex -I- (sind — ricos9)Ey, d = d„ — 9i, and 9i = 
arg(Azi): and y = p/a. 

EDSR coupling DE^^ax in the qubit subspace arises 
from oscillating in-plane electric helds Ey = Ey (x cos dy -|- 
ysindy). We obtain the effective EDSR interactions 
DE\\(jx in the qubit basis using a Schrieffer-Wolff trans¬ 
formation. Expanding in rji = \Xzi\/{si — e„) and rju = 
\Xzu\/{£i — £u) and grouping terms according to powers 
in the electric helds Ei and E 2 , we obtain 


_2\Xzu\Re{X*EiXz2) 

{ei - £„)2 

2Ke{ElXzi + EIXz2) 

-\-a - 

, ^2\Xzu\^e{ElE2) 

+a - 7 --. 


(A24) 


5hu: = -\{^^^)\\Xzi\ - |Az„| cos(2d„ -9i+ d„)). 

(A21) 

The exact solution in (main text, Eig. 2C,D) shows that 
all higher order corrections (including to the ap¬ 

proximate solution presented in the main text do not 
qualitatively modify the qubit frequency. 


The term proportional to E^Xzi + ElXz 2 is linear in the 
electric held and dehnes the electric dipole spin resonance 
term. Substituting Xzi, Xz 2 , Ei and E 2 we obtain the 
EDSR matrix element 

D = a (cos(do — dy) -I- r7sin(do — dy)). (A25) 
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SPIN-DEPENDENT DIPOLE-DIPOLE 
INTERACTION 

Because of the spin-orbit interaction, the electric 
dipole moment in each acceptor couples to spin. As a 
result, two qubits interacting only via mutual Coulomb 
repulsion experience a spin-dependent interaction resem¬ 
bling a magnetic dipole-dipole interaction. Here we de¬ 
termine this interaction in the coupled-qubit subspace 

I-), I—h), |-l—), !++)• For the total Hamiltonian we 

have where is the single¬ 

acceptor Hamiltonian for qubit i = 1,2, and Vi 2 (ri — 
r 2 ) = e^/ 47 re|ri — r 2 | is the electrostatic interaction be¬ 
tween the qubits. 

We work in the tensor 16 x 16 tensor product sub¬ 
space of the two qubits \mn) = (g) |n^) where 

m G { —,-b, M—, u-|-} and n G ,-b, u—, m-I-}, explicitly 
ignoring anti-symmetrization, i.e., assuming spatial over¬ 
laps are negligible. Without the Coulomb interaction, the 
Hamiltonian is 


The combined effect of the Coulomb and Zeeman inter¬ 
actions can be projected into the coupled-qubit subspace 
I ), I—b) , |H—) , I++) using a Schrieffer-Wolff trans¬ 
formation. Working out the effective interaction, to sec¬ 
ond order in off-diagonal terms and zeroth order in \\zu\ 
and \\zi\, gives a spin-independent shift to all levels 
which does not influence the qubit physics. A transfor¬ 
mation to third order in off-diagonal terms in the full 
16 X 16 space finds the spin-spin interaction of interest, 

/ 0 0 0 1 \ 



0 0 10 

Rdd-Jxx 0 10 0 ’ 

\ 1 0 0 0 / 

(A31) 

which is an Ising type spin-spin interaction, 
nience this can be re-written as 

For conve- 

Hdd 

— Jxx^lx^2x 

= Jxx{<^l+ E (7i-){a2-\- + (72-), 

(A32) 

(A33) 


{mn\Hlp + Hlp\m'n') = {m\Hlp\m') 5n 
+ {n\Hlp\n') 5^, 


(A26) 


where is given by Equation (A13). Meanwhile, in 


the direct product subspace the two-qubit Coulomb in¬ 
teraction is 




' , 3e^^L(ri)^n(i‘2)4'm'(ri)4'„'(r2) 


= / drfdr^- 


47re|ri - r 2 | 


(A27) 


When the separation R 12 between the acceptors is large 
compared to dipole moments {Svi) of the system, we may 
use the multi-pole expansion of the Coulomb interaction 


in Equation (A27). The lowest-order non-zero term is 

o 2 


(mn\V 




> = 


47rei?f2 - 


Ri2 • {6r2), 


-3((iJri)w ■ R'i 2 )((^i’ 2 )„m/ - 1 ^ 12 ) 

(A28) 


where 


= J dr^Vi - Ri)4'l(D)4'n'(i’*), (A29) 

Ri is the position of qubit ion i, and R 12 = Ri — 
R 2 . The Coulomb interaction is now a product of 
single-hole dipole matrix elements known in the basis 


{ —,-b, M—, u-b} from Equation (A23). The total Hamil¬ 
tonian is 

{mn\H^\m'n') = {m\H^p\m') Smm' + {n\H^p\n') Snn' 

-b (mn|F^^|m'n') (A30) 


where a± = ajx ± icrjy is the raising/lowering operator 
for qubit j. For R 12 = xi?i 2 , we have contributions Jxx,x 
and Jxx,y from the x and y oriented dipoles, respectively 
such that the total coupling is Jxx = —‘2Jxx,x + Jxx,y, 
where 

T _ 4 (g 2 p|Azi| -b qifj.\Xz 2 \)^ 

qix = asin( 0 ), q 2 x = —acos{9), qiy = acos{0), and 
q 2 y = asin(0). Substituting |Azi| = ^|A^o| cos(—0/2— 
do), |Az 2 | = llAzol sin(- 0/2 - 0 o), we obtain 

^ _ a^|Azopcos2 0o 

Jxx,x — „ r->ri / \o ? 3-nd (A35) 


47rei?3(£^ - ez)2 ’ 

a^lAzopsin^ 0 o 


’^=ox,y 4TreR3(e^ - £ 1 )^’ 


(A36) 


such that 


^ _ /a|Azo^ 2 (- 2 cos 2 0 o-bsin^ 0 o) 

7 T 7 A 3 iA37j 




AireR^ 




(A38) 


PHONON-INDUCED SPIN RELAXATION 

The relaxation from |n') to \n) via emission of a 
phonon with energy huiqs = hvsqa can be determined 
from Fermi’s golden rule. 


1 


T 

rj. 


^ ^ I 5 ^q) | 

xS{En-E'^-hujq,) (A39) 


















where s = are the phonon polarizations, qs is 

the phonon wavevector, H^js = J2i,j Dijeijs is the 
electron-phonon interaction, and riq is the phonon pop¬ 
ulation. The deformation potential matrices Dij are de¬ 
termined from the Bir-Pikus Hamiltonian 


7 . Dijeijs 


=a' {e 


c.p.) 


~ \l)^xxs + c.p.] 

+ {2d' /Vi)[{Jx, Jy}£xys + C.p.] 


(A40) 


where a', b', and d' are Bir-Pikus deformation 
potentials [55] |S7] and the strain djs = ^{dSRis/drj + 
dSRjs/dri) of the phonon polarization s is determined 
by the displacement [75] 


(5Rs = {-i)i 


n 


2NV pujq 


=«s(“qs+»qs)exp(tqs-r), (A41) 


where e^s is the normalized phonon polarization 
vector[74]. N is the number of unit cells, V is the 
unit cell volume, NV = is the crystal volume, p 
is the mass density, and (oqs) creates (destroys) 
a phonon of wavevector qg and polarization s. For 
(n'|i7ij exp(iq • r)|n) we use the dipole approximation 
{n'\Dij{l + ic[-r + .. .)\n) « {n'\Dij\n), which is ap¬ 
propriate since qa ~ 10“^ where q is the the phonon 
wavevector and a ~ 1 nm is the Bohr radius. At low 
temperatures T hui/k ^ 0.7 K we obtain 

{hujf 


1 


20h‘^TTp L 


^1 {n\D,,\n) 


+ 


ZvV 


-f^|(n'|A,|n)p(^ + ^)] (A42) 

<■ * 


where vi and Vt are the longitudinal and transverse sound 
velocities. For our qubit, the relaxation from |-|-) to |—) 
can be evaluated using a Schrieffer-Wolff transformation 
of the non-diagonal elements oi H = Hop + Heijs, where 
Hsijs = Deijs^ijaWto- We determine the cou¬ 

pling to lowest order in Hzo and H^ijs while treating 
He, ion, Hif and He exactly within the 4x4 subspace. 
We obtain 


1 {fiujy / Huj N2r9 , 2 /- 
fi ~ 20h^TTp\4pEj ^ 32 ° 


SviJ 




16 


\3vf 


(A43) 


at the sweet spot. For comparison, we also evaluated the 
phonon-mediated transition rate from the first excited 
level | —^) to the ground state | —|) in bulk unstrained 
silicon, which is allowed even in zero magnetic field, since 
| —^) and | —I) are not time-reversal symmetric. For Bq 
along [001] directions we obtain 


1 




■a 


3uf 


2d‘ 


r/2 


20p'Khd \Vf 

in agreement with the low-temperature result in ref. 


(A44) 


DEPHASING FROM ELECTRIC FIELD 
FLUCTUATIONS 

Dephasing of a spin-orbit qubit occurs due to random 
fluctuations hSuj(t) in the qubit level energy splitting. In 
this section we outline expected dephasing rates associ¬ 
ated with realistic parameters for the qubit. 


Charge noise from trap charging/discharging 


First, we estimate the dephasing due to a single fluc¬ 
tuating charge trap (random telegraph signal), assum¬ 
ing that a charge trap can only fluctuate when it is in 
tunneling proximity to carrier reservoir or gate. The 
dephasing rate is given by (T^)”^ = {6oj)'^ts/2, where 
Shoj is the qubit energy shift when the trap is charged, 
and Ts is the switching time [36). Together with its im¬ 
age in the gate, a dipole e£d is created, resulting in a 
dipole potential Vd = eidzj^TT^x^ + when the de¬ 

fect is a distance R = icx + zz away from the accep¬ 
tor qubit. The electric field for the dipole potential is 
(5E = -VVd = ±SEs, + z5E^. 

The change Stiuj of the Larmor frequency is readily cal¬ 


culated from the extended qubit Hamiltonian (Eq. A231, 


where 5Ex,y couples off-diagonally through the aEi + Xzi 
terms, and SE^ on-diagonally via the explicit dependence 
of A{E^) and \XziiEz)\ on E^. Equivalently, the 2x2 
qubit model can be used, and it is easy to show that 


2D^6E^ „ 

Suj « uj{Ez + SEz) — uj{Ez) + —. „ ( • 

h^w{Ez) 


(A45) 


For the estimate of dephasing in the main text we as¬ 
sumed a: = 50 nm, and z = 5 nm giving 5Ex ~ 600 V/m 
and 5Ez « 2000 V/m. 


Gate electric field noise 


Gate electric field noise is modeled as a white John¬ 
son voltage noise u„(t) process applied across the gates 
producing a randomly fluctuating field of order E{t) « 
Vnit)/dg where dg « 20 nm is the shortest envisioned 
distance between gates. The dephasing rate is given by 
(T|)“^ = S{uj = 0) where S{oj) is the power spectral den¬ 
sity of the random process 6u;{t). Further approximation 
to Slu from above gives 


Suj{t) 


duj{Ez) Vnjt) 20"^ / u„(f) \2 

dEz dg h?io \ dg ) 


(A46) 


We find that the first term is much larger for our qubit 
and dominates the power spectral density of duj{t), 


S{u}) « ( 


dcjjEz) 1 
dEz dg 



(A47) 
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The quantity duj{Ez)/dEz is determined from our ana¬ 
lytic model. Substituting S^iuj) = AksTR as the white 
noise power spectrum, T = 1 K and R = 50 ohms, we 
find that the dephasing rate due to intrinsic John¬ 

son voltage noise on the gate is 10 ^ times smaller than 
the dephasing rate from charge noise. 

NUMERICAL KOHN LUTTINGER 
CALCULATIONS 

All theory preditions in the main text were compared 
with a numerical solution of the acceptor Hamiltonian in¬ 
cluding the full spatial dependence of J/ (k) . The Hamil¬ 
tonian H = JJlk + He, ion, where JJlk = Hik) + + 

JJion(r) -I- Hii{z) -\- He + Hz, was computed in the 6 x 6 
representation of valence band Bloch states |T, mj)[45], 
\fpmj) = \J,mj), where Fj,mj{r) are en¬ 

velope functions. The first step is to numerically diago- 


nalize Hek 

using H{k) = 





^ P E Q 

-S 

R 

0 

- 1 s 

V2R 

\ 

-5t 

P-Q 

0 

s 

-y/2R 




0 

P-Q 

s 

\A-5' 

y/2Q 


0 

R) 


PEQ 

-V2R) 




-V2Q 


-V2R 

P E Aso 

0 


\ V2R 

/ 3 'ct 
V 2*^ 


- 1 S' 

0 

P E Aso 

/ 




(A48j 

in the |,i 

\ 13 1 

/’ 12 ’ 2 

'I P 

/’ 1 2 ’ 

1\ 13 

2/’ 1 2 ’ 

3\ 111 

2/’ 12 ’ 2 

'I P ^ 

/’ 12 ’ 2 



basis, with 


-P — (A49) 

Q=^^l2{kl + kl-2kl), (A50) 

R = ~ ^ 2i-i3kz:ky), (A51) 

h2 

S = - - 2V3^3(kx - iky)kz, (A52) 

ZTTIq 

^^ion(i') = -e^/47rer. He = -eE^z, H^f = Uo&{-z) with 
Uq ^ oo (using a Dirichlet boundary condition), and Hz 
and TJe from the main text. Here, kj, = i-§^, ky = 
kz = i-§^, and tij are the strains. The dielectric constant 
e = 11.4eo, LK parameters 71 = 4.28, 72 = 0.375 and 
73 = 1.44, and Aso = 44 meV, Bir-Pikus deformation 
potentials b' = —1.42 eV and d! = —3.7 eV, and Lande g- 
factors gi = 1.07 and g 2 = 0.033 were obtained from the 
literature [m |M1 ESI llH]. The hydrostatic deformation 
potential a' is not needed since it acts on the valence 
bands with an identity operator. 

The Hamiltonian Hlk was diagonalized by finite differ¬ 
ences in Cartesian coordinates for different fields Ez and 


acceptor depths d, to obtain the eigen-energies (eigen¬ 
states) Elk (I'I'lk)) of 77 lk- The inversion asymmetry- 
induced mixing of states outside the 4x4 subspace is 
explicitly included due to the explicit representation of 
electric potentials of the ion, interface, and applied gate 
field, in the Hamiltonian. The effect of the Td symmetry 
terms is absent in I^lk) because the local field of the ion 
is only present in the unit cell containing the ion, which 
lies outside the validity of the LK approach [551157] . 

The Td symmetry term was included in the spectrum 
and eigenstates by diagonalizing H^k E He, ion in the ba¬ 
sis I^lk)- The p values in HE,ion were taken to repro¬ 
duce LH-HH couplings known from experiments [58] for 
calculations on bulk acceptors. We reduced p by a mul¬ 
tiplicative prefactor due to the interface and gate field, 
which was determined from the reduction of the charge 
in density at the ion in I'I'lk)- 

We diagonalized H^k + HE,ion using the first 8 states 
from JJlk- We note that states 5 and 6 contribute 
very little; even though these predominantly split-off hole 
states have a non-zero density at the ion (non-zero TJ 
symmetry correction), they are already 20 meV away 
(compared to the ~ 0.2 meV HH-LH level splitting). The 
even higher energy states 7 and higher may be neglected 
a priori at low fields, because they have even higher ener¬ 
gies, but more importantly, since their principal quantum 
number is 2 or higher. Consequently, the density at the 
ion, and hence the Td symmetry coupling, is even smaller. 

Hif + He'- symmetry and parameters for analytic 
model 

The states IJ^lk) experience inversion asymmetry due 
to explicit presence of the interface and gate field. We 
have verified that inversion asymmetry-induced mixing 
of excited states into the 4x4 subspace does not mod¬ 
ify the symmetry of He from the form in the main 
text, for different acceptor depths, by evaluating He = 
(4'm'j,LK|eJE • using eigenstates |Tlk) calcu¬ 

lated non-perturbatively, that is, with applied potentials 
including the ion, interface, and gate electric field. This 
was checked for all static fields Ez in the main text, and 
for fields JE for all three orthogonal directions. Values 
for Aif J- A{Ez) were obtained from eigen-energies, while 
values of a were obtained from the off-diagonal elements 
coupling the LH and HH blocks. 

Larmor frequency and EDSR: numerical results 

The LH-HH mixing caused by p in He, ion changes 
the Larmor frequency hw of our spin qubit |±) dra¬ 
matically compared with |4 ^±lk)- Numerical calcula¬ 
tions of hui (Fig. 2C,D, blue squares) are in agreement 
with the approximate (Fig. 2C,D, black lines) and an- 
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alytic solutions (Fig. 2C,D, green lines) to the 4x4 
model. In particular, our numerical calculations show 
that mixing of states outside the 4x4 subspace, and 
the inclusion of cubic Lande g-factor g 2 with values ob¬ 
tained in experiments, does not qualitatively change the 
behaviour of the qubit splitting. Therefore, the main 
electric-field dependence on electric fields is dominated 
by the LH-HH mixing inside the 4x4 subspace. The 
numerical results for EDSR (Fig. 2E,F, orange squares) 
were obtained by transforming our numerical represen¬ 
tations of He and He, ion into the basis of H, using 
(-|-|eE • r|—) = U'^{He + HE,ion)U\ where U' is the uni¬ 
tary matrix that diagonalizes Hek + He, ion in the basis 
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